R Notebook
This is an R Markdown Notebook. When you execute code within the notebook, the results appear beneath the code.
Try executing this chunk by clicking the Run button within the chunk or by placing your cursor inside it and pressing Ctrl+Shift+Enter.
source("./tianfengRwrappers.R")# plan("multiprocess",workers = 8)读取结果
human_coronary <- readRDS("human_coronary.rds")
CA_dataset1 <- readRDS("CA_dataset1.rds")
CA_dataset2 <- readRDS("CA_dataset2.rds") #已经经过分组处理了修改分群
umapplot(CA_dataset2, split.by = "sample")
umapplot(CA_dataset2,group.by = "groups", split.by = "conditions")
multi_featureplot(c("HEY1","GJA5","SEMA3G","CXCL12","SOX17","CDH5","PECAM1"),CA_dataset2)
multi_featureplot(c("ACKR1","PLVAP","ITGA6","PECAM1"),CA_dataset2)
# table(CA_dataset2$sample)
# 关注cluster 3-6 EC
# cluster 7 基质细胞
multi_featureplot(c("LYZ","PTPRC","CD69","EPCAM","CDH1","PDGFRB","COL1A2","PECAM1","CLDN5"),CA_dataset2)
Dotplot(c("LUM","MMP2","MGP","DCN","MYH11","ACTA2","CNN1","TAGLN"),CA_dataset2) #cluster 7 细胞可以被认为是modulated SMCs
CA_dataset2 <- AddModuleScore(CA_dataset2,list(c("LUM","MMP2","MGP","DCN")))
CA_dataset2 <- AddModuleScore(CA_dataset2,list(c("MYH11","ACTA2","CNN1","TAGLN")))
multi_featureplot(c("Cluster1","LUM","ACTA2","TAGLN"),CA_dataset2)
multi_featureplot(c("MMP2","GJA4","PECAM1","ACKR1"), CA_dataset2)
multi_featureplot(c("MMP2","GJA4","PECAM1","ACKR1"), human_coronary)
细胞分类
multi_featureplot(c("LYZ","PTPRC","CD69","PDGFRB","COL1A2","PECAM1","CLDN5"), CA_dataset1)
f("LUM",CA_dataset1)
umapplot(CA_dataset1)CA_dataset1 <- CA_dataset1 %>% FindNeighbors(dims = 1:20) %>% FindClusters(resolution = 0.5)Computing nearest neighbor graph
Computing SNNModularity Optimizer version 1.3.0 by Ludo Waltman and Nees Jan van Eck
Number of nodes: 8709
Number of edges: 279176
Running Louvain algorithm...0% 10 20 30 40 50 60 70 80 90 100%
[----|----|----|----|----|----|----|----|----|----|
**************************************************|Maximum modularity in 10 random starts: 0.9035
Number of communities: 13
Elapsed time: 0 seconds
CA_dataset2_markers <- FindAllMarkers(CA_dataset2, logfc.threshold = 0.7, min.diff.pct = 0.3,only.pos = T)
Calculating cluster T cell
| | 0 % ~calculating
|+++ | 4 % ~01s
|+++++ | 9 % ~01s
|+++++++ | 13% ~01s
|+++++++++ | 17% ~01s
|+++++++++++ | 22% ~01s
|++++++++++++++ | 26% ~01s
|++++++++++++++++ | 30% ~01s
|++++++++++++++++++ | 35% ~01s
|++++++++++++++++++++ | 39% ~00s
|++++++++++++++++++++++ | 43% ~00s
|++++++++++++++++++++++++ | 48% ~00s
|+++++++++++++++++++++++++++ | 52% ~00s
|+++++++++++++++++++++++++++++ | 57% ~00s
|+++++++++++++++++++++++++++++++ | 61% ~00s
|+++++++++++++++++++++++++++++++++ | 65% ~00s
|+++++++++++++++++++++++++++++++++++ | 70% ~00s
|+++++++++++++++++++++++++++++++++++++ | 74% ~00s
|++++++++++++++++++++++++++++++++++++++++ | 78% ~00s
|++++++++++++++++++++++++++++++++++++++++++ | 83% ~00s
|++++++++++++++++++++++++++++++++++++++++++++ | 87% ~00s
|++++++++++++++++++++++++++++++++++++++++++++++ | 91% ~00s
|++++++++++++++++++++++++++++++++++++++++++++++++ | 96% ~00s
|++++++++++++++++++++++++++++++++++++++++++++++++++| 100% elapsed=01s Calculating cluster NK cell
| | 0 % ~calculating
|+++ | 5 % ~01s
|+++++ | 10% ~01s
|++++++++ | 15% ~01s
|++++++++++ | 20% ~01s
|+++++++++++++ | 25% ~01s
|+++++++++++++++ | 30% ~00s
|++++++++++++++++++ | 35% ~00s
|++++++++++++++++++++ | 40% ~00s
|+++++++++++++++++++++++ | 45% ~00s
|+++++++++++++++++++++++++ | 50% ~00s
|++++++++++++++++++++++++++++ | 55% ~00s
|++++++++++++++++++++++++++++++ | 60% ~00s
|+++++++++++++++++++++++++++++++++ | 65% ~00s
|+++++++++++++++++++++++++++++++++++ | 70% ~00s
|++++++++++++++++++++++++++++++++++++++ | 75% ~00s
|++++++++++++++++++++++++++++++++++++++++ | 80% ~00s
|+++++++++++++++++++++++++++++++++++++++++++ | 85% ~00s
|+++++++++++++++++++++++++++++++++++++++++++++ | 90% ~00s
|++++++++++++++++++++++++++++++++++++++++++++++++ | 95% ~00s
|++++++++++++++++++++++++++++++++++++++++++++++++++| 100% elapsed=01s Calculating cluster Neural progenitor
| | 0 % ~calculating
|+ | 2 % ~02s
|++ | 3 % ~02s
|+++ | 5 % ~02s
|++++ | 7 % ~02s
|+++++ | 8 % ~02s
|++++++ | 10% ~02s
|++++++ | 12% ~02s
|+++++++ | 14% ~02s
|++++++++ | 15% ~02s
|+++++++++ | 17% ~02s
|++++++++++ | 19% ~02s
|+++++++++++ | 20% ~02s
|++++++++++++ | 22% ~02s
|++++++++++++ | 24% ~02s
|+++++++++++++ | 25% ~01s
|++++++++++++++ | 27% ~01s
|+++++++++++++++ | 29% ~01s
|++++++++++++++++ | 31% ~01s
|+++++++++++++++++ | 32% ~01s
|+++++++++++++++++ | 34% ~01s
|++++++++++++++++++ | 36% ~01s
|+++++++++++++++++++ | 37% ~01s
|++++++++++++++++++++ | 39% ~01s
|+++++++++++++++++++++ | 41% ~01s
|++++++++++++++++++++++ | 42% ~01s
|+++++++++++++++++++++++ | 44% ~01s
|+++++++++++++++++++++++ | 46% ~01s
|++++++++++++++++++++++++ | 47% ~01s
|+++++++++++++++++++++++++ | 49% ~01s
|++++++++++++++++++++++++++ | 51% ~01s
|+++++++++++++++++++++++++++ | 53% ~01s
|++++++++++++++++++++++++++++ | 54% ~01s
|++++++++++++++++++++++++++++ | 56% ~01s
|+++++++++++++++++++++++++++++ | 58% ~01s
|++++++++++++++++++++++++++++++ | 59% ~01s
|+++++++++++++++++++++++++++++++ | 61% ~01s
|++++++++++++++++++++++++++++++++ | 63% ~01s
|+++++++++++++++++++++++++++++++++ | 64% ~01s
|++++++++++++++++++++++++++++++++++ | 66% ~01s
|++++++++++++++++++++++++++++++++++ | 68% ~01s
|+++++++++++++++++++++++++++++++++++ | 69% ~01s
|++++++++++++++++++++++++++++++++++++ | 71% ~01s
|+++++++++++++++++++++++++++++++++++++ | 73% ~01s
|++++++++++++++++++++++++++++++++++++++ | 75% ~01s
|+++++++++++++++++++++++++++++++++++++++ | 76% ~00s
|+++++++++++++++++++++++++++++++++++++++ | 78% ~00s
|++++++++++++++++++++++++++++++++++++++++ | 80% ~00s
|+++++++++++++++++++++++++++++++++++++++++ | 81% ~00s
|++++++++++++++++++++++++++++++++++++++++++ | 83% ~00s
|+++++++++++++++++++++++++++++++++++++++++++ | 85% ~00s
|++++++++++++++++++++++++++++++++++++++++++++ | 86% ~00s
|+++++++++++++++++++++++++++++++++++++++++++++ | 88% ~00s
|+++++++++++++++++++++++++++++++++++++++++++++ | 90% ~00s
|++++++++++++++++++++++++++++++++++++++++++++++ | 92% ~00s
|+++++++++++++++++++++++++++++++++++++++++++++++ | 93% ~00s
|++++++++++++++++++++++++++++++++++++++++++++++++ | 95% ~00s
|+++++++++++++++++++++++++++++++++++++++++++++++++ | 97% ~00s
|++++++++++++++++++++++++++++++++++++++++++++++++++| 98% ~00s
|++++++++++++++++++++++++++++++++++++++++++++++++++| 100% elapsed=02s Calculating cluster Stromal cell
| | 0 % ~calculating
|+ | 2 % ~04s
|++ | 4 % ~04s
|+++ | 5 % ~04s
|++++ | 7 % ~04s
|+++++ | 9 % ~04s
|++++++ | 11% ~04s
|+++++++ | 12% ~03s
|++++++++ | 14% ~03s
|+++++++++ | 16% ~03s
|+++++++++ | 18% ~04s
|++++++++++ | 20% ~04s
|+++++++++++ | 21% ~04s
|++++++++++++ | 23% ~04s
|+++++++++++++ | 25% ~04s
|++++++++++++++ | 27% ~04s
|+++++++++++++++ | 29% ~03s
|++++++++++++++++ | 30% ~03s
|+++++++++++++++++ | 32% ~03s
|+++++++++++++++++ | 34% ~03s
|++++++++++++++++++ | 36% ~03s
|+++++++++++++++++++ | 38% ~03s
|++++++++++++++++++++ | 39% ~03s
|+++++++++++++++++++++ | 41% ~03s
|++++++++++++++++++++++ | 43% ~03s
|+++++++++++++++++++++++ | 45% ~03s
|++++++++++++++++++++++++ | 46% ~03s
|+++++++++++++++++++++++++ | 48% ~03s
|+++++++++++++++++++++++++ | 50% ~03s
|++++++++++++++++++++++++++ | 52% ~02s
|+++++++++++++++++++++++++++ | 54% ~02s
|++++++++++++++++++++++++++++ | 55% ~02s
|+++++++++++++++++++++++++++++ | 57% ~02s
|++++++++++++++++++++++++++++++ | 59% ~02s
|+++++++++++++++++++++++++++++++ | 61% ~02s
|++++++++++++++++++++++++++++++++ | 62% ~02s
|+++++++++++++++++++++++++++++++++ | 64% ~02s
|++++++++++++++++++++++++++++++++++ | 66% ~02s
|++++++++++++++++++++++++++++++++++ | 68% ~02s
|+++++++++++++++++++++++++++++++++++ | 70% ~01s
|++++++++++++++++++++++++++++++++++++ | 71% ~01s
|+++++++++++++++++++++++++++++++++++++ | 73% ~01s
|++++++++++++++++++++++++++++++++++++++ | 75% ~01s
|+++++++++++++++++++++++++++++++++++++++ | 77% ~01s
|++++++++++++++++++++++++++++++++++++++++ | 79% ~01s
|+++++++++++++++++++++++++++++++++++++++++ | 80% ~01s
|++++++++++++++++++++++++++++++++++++++++++ | 82% ~01s
|++++++++++++++++++++++++++++++++++++++++++ | 84% ~01s
|+++++++++++++++++++++++++++++++++++++++++++ | 86% ~01s
|++++++++++++++++++++++++++++++++++++++++++++ | 88% ~01s
|+++++++++++++++++++++++++++++++++++++++++++++ | 89% ~00s
|++++++++++++++++++++++++++++++++++++++++++++++ | 91% ~00s
|+++++++++++++++++++++++++++++++++++++++++++++++ | 93% ~00s
|++++++++++++++++++++++++++++++++++++++++++++++++ | 95% ~00s
|+++++++++++++++++++++++++++++++++++++++++++++++++ | 96% ~00s
|++++++++++++++++++++++++++++++++++++++++++++++++++| 98% ~00s
|++++++++++++++++++++++++++++++++++++++++++++++++++| 100% elapsed=05s Calculating cluster EC
| | 0 % ~calculating
|+ | 1 % ~02s
|++ | 3 % ~02s
|+++ | 4 % ~02s
|+++ | 6 % ~02s
|++++ | 7 % ~02s
|+++++ | 8 % ~02s
|+++++ | 10% ~02s
|++++++ | 11% ~02s
|+++++++ | 13% ~02s
|++++++++ | 14% ~02s
|++++++++ | 15% ~02s
|+++++++++ | 17% ~02s
|++++++++++ | 18% ~02s
|++++++++++ | 20% ~02s
|+++++++++++ | 21% ~02s
|++++++++++++ | 23% ~02s
|++++++++++++ | 24% ~02s
|+++++++++++++ | 25% ~02s
|++++++++++++++ | 27% ~02s
|+++++++++++++++ | 28% ~02s
|+++++++++++++++ | 30% ~02s
|++++++++++++++++ | 31% ~02s
|+++++++++++++++++ | 32% ~02s
|+++++++++++++++++ | 34% ~02s
|++++++++++++++++++ | 35% ~02s
|+++++++++++++++++++ | 37% ~02s
|++++++++++++++++++++ | 38% ~02s
|++++++++++++++++++++ | 39% ~02s
|+++++++++++++++++++++ | 41% ~02s
|++++++++++++++++++++++ | 42% ~02s
|++++++++++++++++++++++ | 44% ~01s
|+++++++++++++++++++++++ | 45% ~01s
|++++++++++++++++++++++++ | 46% ~01s
|++++++++++++++++++++++++ | 48% ~01s
|+++++++++++++++++++++++++ | 49% ~01s
|++++++++++++++++++++++++++ | 51% ~01s
|+++++++++++++++++++++++++++ | 52% ~01s
|+++++++++++++++++++++++++++ | 54% ~01s
|++++++++++++++++++++++++++++ | 55% ~01s
|+++++++++++++++++++++++++++++ | 56% ~01s
|+++++++++++++++++++++++++++++ | 58% ~01s
|++++++++++++++++++++++++++++++ | 59% ~01s
|+++++++++++++++++++++++++++++++ | 61% ~01s
|+++++++++++++++++++++++++++++++ | 62% ~01s
|++++++++++++++++++++++++++++++++ | 63% ~01s
|+++++++++++++++++++++++++++++++++ | 65% ~01s
|++++++++++++++++++++++++++++++++++ | 66% ~01s
|++++++++++++++++++++++++++++++++++ | 68% ~01s
|+++++++++++++++++++++++++++++++++++ | 69% ~01s
|++++++++++++++++++++++++++++++++++++ | 70% ~01s
|++++++++++++++++++++++++++++++++++++ | 72% ~01s
|+++++++++++++++++++++++++++++++++++++ | 73% ~01s
|++++++++++++++++++++++++++++++++++++++ | 75% ~01s
|+++++++++++++++++++++++++++++++++++++++ | 76% ~01s
|+++++++++++++++++++++++++++++++++++++++ | 77% ~01s
|++++++++++++++++++++++++++++++++++++++++ | 79% ~01s
|+++++++++++++++++++++++++++++++++++++++++ | 80% ~01s
|+++++++++++++++++++++++++++++++++++++++++ | 82% ~00s
|++++++++++++++++++++++++++++++++++++++++++ | 83% ~00s
|+++++++++++++++++++++++++++++++++++++++++++ | 85% ~00s
|+++++++++++++++++++++++++++++++++++++++++++ | 86% ~00s
|++++++++++++++++++++++++++++++++++++++++++++ | 87% ~00s
|+++++++++++++++++++++++++++++++++++++++++++++ | 89% ~00s
|++++++++++++++++++++++++++++++++++++++++++++++ | 90% ~00s
|++++++++++++++++++++++++++++++++++++++++++++++ | 92% ~00s
|+++++++++++++++++++++++++++++++++++++++++++++++ | 93% ~00s
|++++++++++++++++++++++++++++++++++++++++++++++++ | 94% ~00s
|++++++++++++++++++++++++++++++++++++++++++++++++ | 96% ~00s
|+++++++++++++++++++++++++++++++++++++++++++++++++ | 97% ~00s
|++++++++++++++++++++++++++++++++++++++++++++++++++| 99% ~00s
|++++++++++++++++++++++++++++++++++++++++++++++++++| 100% elapsed=03s Calculating cluster Macrophage
| | 0 % ~calculating
|+ | 2 % ~04s
|++ | 3 % ~04s
|+++ | 5 % ~04s
|++++ | 7 % ~04s
|+++++ | 8 % ~04s
|++++++ | 10% ~04s
|++++++ | 12% ~04s
|+++++++ | 14% ~04s
|++++++++ | 15% ~04s
|+++++++++ | 17% ~04s
|++++++++++ | 19% ~04s
|+++++++++++ | 20% ~03s
|++++++++++++ | 22% ~03s
|++++++++++++ | 24% ~03s
|+++++++++++++ | 25% ~03s
|++++++++++++++ | 27% ~03s
|+++++++++++++++ | 29% ~03s
|++++++++++++++++ | 31% ~03s
|+++++++++++++++++ | 32% ~03s
|+++++++++++++++++ | 34% ~03s
|++++++++++++++++++ | 36% ~03s
|+++++++++++++++++++ | 37% ~03s
|++++++++++++++++++++ | 39% ~03s
|+++++++++++++++++++++ | 41% ~03s
|++++++++++++++++++++++ | 42% ~03s
|+++++++++++++++++++++++ | 44% ~02s
|+++++++++++++++++++++++ | 46% ~02s
|++++++++++++++++++++++++ | 47% ~02s
|+++++++++++++++++++++++++ | 49% ~02s
|++++++++++++++++++++++++++ | 51% ~02s
|+++++++++++++++++++++++++++ | 53% ~02s
|++++++++++++++++++++++++++++ | 54% ~02s
|++++++++++++++++++++++++++++ | 56% ~02s
|+++++++++++++++++++++++++++++ | 58% ~02s
|++++++++++++++++++++++++++++++ | 59% ~02s
|+++++++++++++++++++++++++++++++ | 61% ~02s
|++++++++++++++++++++++++++++++++ | 63% ~02s
|+++++++++++++++++++++++++++++++++ | 64% ~02s
|++++++++++++++++++++++++++++++++++ | 66% ~02s
|++++++++++++++++++++++++++++++++++ | 68% ~01s
|+++++++++++++++++++++++++++++++++++ | 69% ~01s
|++++++++++++++++++++++++++++++++++++ | 71% ~01s
|+++++++++++++++++++++++++++++++++++++ | 73% ~01s
|++++++++++++++++++++++++++++++++++++++ | 75% ~01s
|+++++++++++++++++++++++++++++++++++++++ | 76% ~01s
|+++++++++++++++++++++++++++++++++++++++ | 78% ~01s
|++++++++++++++++++++++++++++++++++++++++ | 80% ~01s
|+++++++++++++++++++++++++++++++++++++++++ | 81% ~01s
|++++++++++++++++++++++++++++++++++++++++++ | 83% ~01s
|+++++++++++++++++++++++++++++++++++++++++++ | 85% ~01s
|++++++++++++++++++++++++++++++++++++++++++++ | 86% ~01s
|+++++++++++++++++++++++++++++++++++++++++++++ | 88% ~01s
|+++++++++++++++++++++++++++++++++++++++++++++ | 90% ~00s
|++++++++++++++++++++++++++++++++++++++++++++++ | 92% ~00s
|+++++++++++++++++++++++++++++++++++++++++++++++ | 93% ~00s
|++++++++++++++++++++++++++++++++++++++++++++++++ | 95% ~00s
|+++++++++++++++++++++++++++++++++++++++++++++++++ | 97% ~00s
|++++++++++++++++++++++++++++++++++++++++++++++++++| 98% ~00s
|++++++++++++++++++++++++++++++++++++++++++++++++++| 100% elapsed=04s Calculating cluster B cell
| | 0 % ~calculating
|+++++ | 9 % ~00s
|++++++++++ | 18% ~00s
|++++++++++++++ | 27% ~00s
|+++++++++++++++++++ | 36% ~00s
|+++++++++++++++++++++++ | 45% ~00s
|++++++++++++++++++++++++++++ | 55% ~00s
|++++++++++++++++++++++++++++++++ | 64% ~00s
|+++++++++++++++++++++++++++++++++++++ | 73% ~00s
|+++++++++++++++++++++++++++++++++++++++++ | 82% ~00s
|++++++++++++++++++++++++++++++++++++++++++++++ | 91% ~00s
|++++++++++++++++++++++++++++++++++++++++++++++++++| 100% elapsed=00s Calculating cluster Mast cell
| | 0 % ~calculating
|++ | 2 % ~01s
|+++ | 5 % ~01s
|++++ | 7 % ~01s
|+++++ | 10% ~01s
|++++++ | 12% ~01s
|++++++++ | 14% ~01s
|+++++++++ | 17% ~01s
|++++++++++ | 19% ~01s
|+++++++++++ | 21% ~01s
|++++++++++++ | 24% ~01s
|++++++++++++++ | 26% ~01s
|+++++++++++++++ | 29% ~01s
|++++++++++++++++ | 31% ~01s
|+++++++++++++++++ | 33% ~01s
|++++++++++++++++++ | 36% ~01s
|++++++++++++++++++++ | 38% ~01s
|+++++++++++++++++++++ | 40% ~01s
|++++++++++++++++++++++ | 43% ~01s
|+++++++++++++++++++++++ | 45% ~01s
|++++++++++++++++++++++++ | 48% ~01s
|+++++++++++++++++++++++++ | 50% ~01s
|+++++++++++++++++++++++++++ | 52% ~01s
|++++++++++++++++++++++++++++ | 55% ~01s
|+++++++++++++++++++++++++++++ | 57% ~01s
|++++++++++++++++++++++++++++++ | 60% ~01s
|+++++++++++++++++++++++++++++++ | 62% ~01s
|+++++++++++++++++++++++++++++++++ | 64% ~01s
|++++++++++++++++++++++++++++++++++ | 67% ~00s
|+++++++++++++++++++++++++++++++++++ | 69% ~00s
|++++++++++++++++++++++++++++++++++++ | 71% ~00s
|+++++++++++++++++++++++++++++++++++++ | 74% ~00s
|+++++++++++++++++++++++++++++++++++++++ | 76% ~00s
|++++++++++++++++++++++++++++++++++++++++ | 79% ~00s
|+++++++++++++++++++++++++++++++++++++++++ | 81% ~00s
|++++++++++++++++++++++++++++++++++++++++++ | 83% ~00s
|+++++++++++++++++++++++++++++++++++++++++++ | 86% ~00s
|+++++++++++++++++++++++++++++++++++++++++++++ | 88% ~00s
|++++++++++++++++++++++++++++++++++++++++++++++ | 90% ~00s
|+++++++++++++++++++++++++++++++++++++++++++++++ | 93% ~00s
|++++++++++++++++++++++++++++++++++++++++++++++++ | 95% ~00s
|+++++++++++++++++++++++++++++++++++++++++++++++++ | 98% ~00s
|++++++++++++++++++++++++++++++++++++++++++++++++++| 100% elapsed=01s 基质细胞
ds0 <- readRDS("ds0.rds")
ds1 <- readRDS("ds1.rds")
ds2 <- readRDS("ds2.rds")去除边角的离群细胞
select.cells <- CellSelector(plot = DimPlot(ds0, reduction = "umap"))
Listening on http://127.0.0.1:3689
ds0 <- subset(ds0, cell = select.cells)
umapplot(ds0)
细分亚群
ds2_PA_markers <- FindAllMarkers(ds2_PA, logfc.threshold = 0.7, min.diff.pct = 0.3,only.pos = T)
Calculating cluster Fibromyocyte
| | 0 % ~calculating
|+++ | 6 % ~00s
|++++++ | 12% ~00s
|+++++++++ | 18% ~00s
|++++++++++++ | 24% ~00s
|+++++++++++++++ | 29% ~00s
|++++++++++++++++++ | 35% ~00s
|+++++++++++++++++++++ | 41% ~00s
|++++++++++++++++++++++++ | 47% ~00s
|+++++++++++++++++++++++++++ | 53% ~00s
|++++++++++++++++++++++++++++++ | 59% ~00s
|+++++++++++++++++++++++++++++++++ | 65% ~00s
|++++++++++++++++++++++++++++++++++++ | 71% ~00s
|+++++++++++++++++++++++++++++++++++++++ | 76% ~00s
|++++++++++++++++++++++++++++++++++++++++++ | 82% ~00s
|+++++++++++++++++++++++++++++++++++++++++++++ | 88% ~00s
|++++++++++++++++++++++++++++++++++++++++++++++++ | 94% ~00s
|++++++++++++++++++++++++++++++++++++++++++++++++++| 100% elapsed=00s Calculating cluster SMC1
| | 0 % ~calculating
|++++ | 7 % ~00s
|++++++++ | 14% ~00s
|+++++++++++ | 21% ~00s
|+++++++++++++++ | 29% ~00s
|++++++++++++++++++ | 36% ~00s
|++++++++++++++++++++++ | 43% ~00s
|+++++++++++++++++++++++++ | 50% ~00s
|+++++++++++++++++++++++++++++ | 57% ~00s
|+++++++++++++++++++++++++++++++++ | 64% ~00s
|++++++++++++++++++++++++++++++++++++ | 71% ~00s
|++++++++++++++++++++++++++++++++++++++++ | 79% ~00s
|+++++++++++++++++++++++++++++++++++++++++++ | 86% ~00s
|+++++++++++++++++++++++++++++++++++++++++++++++ | 93% ~00s
|++++++++++++++++++++++++++++++++++++++++++++++++++| 100% elapsed=00s Calculating cluster SMC2
| | 0 % ~calculating
|+++++ | 9 % ~00s
|++++++++++ | 18% ~00s
|++++++++++++++ | 27% ~00s
|+++++++++++++++++++ | 36% ~00s
|+++++++++++++++++++++++ | 45% ~00s
|++++++++++++++++++++++++++++ | 55% ~00s
|++++++++++++++++++++++++++++++++ | 64% ~00s
|+++++++++++++++++++++++++++++++++++++ | 73% ~00s
|+++++++++++++++++++++++++++++++++++++++++ | 82% ~00s
|++++++++++++++++++++++++++++++++++++++++++++++ | 91% ~00s
|++++++++++++++++++++++++++++++++++++++++++++++++++| 100% elapsed=00s 
write.csv(ds2_AC_markers,"ds2_AC_markers.csv")
write.csv(ds2_PA_markers,"ds2_PA_markers.csv")
LS0tCnRpdGxlOiAiUiBOb3RlYm9vayIKb3V0cHV0OiBodG1sX25vdGVib29rCi0tLQoKVGhpcyBpcyBhbiBbUiBNYXJrZG93bl0oaHR0cDovL3JtYXJrZG93bi5yc3R1ZGlvLmNvbSkgTm90ZWJvb2suIFdoZW4geW91IGV4ZWN1dGUgY29kZSB3aXRoaW4gdGhlIG5vdGVib29rLCB0aGUgcmVzdWx0cyBhcHBlYXIgYmVuZWF0aCB0aGUgY29kZS4gCgpUcnkgZXhlY3V0aW5nIHRoaXMgY2h1bmsgYnkgY2xpY2tpbmcgdGhlICpSdW4qIGJ1dHRvbiB3aXRoaW4gdGhlIGNodW5rIG9yIGJ5IHBsYWNpbmcgeW91ciBjdXJzb3IgaW5zaWRlIGl0IGFuZCBwcmVzc2luZyAqQ3RybCtTaGlmdCtFbnRlciouIAoKYGBge3J9CnNvdXJjZSgiLi90aWFuZmVuZ1J3cmFwcGVycy5SIikKIyBwbGFuKCJtdWx0aXByb2Nlc3MiLHdvcmtlcnMgPSA4KQpgYGAKCiMg6K+75Y+W57uT5p6cCmBgYHtyfQpodW1hbl9jb3JvbmFyeSA8LSByZWFkUkRTKCJodW1hbl9jb3JvbmFyeS5yZHMiKQpDQV9kYXRhc2V0MSA8LSByZWFkUkRTKCJDQV9kYXRhc2V0MS5yZHMiKQpDQV9kYXRhc2V0MiA8LSByZWFkUkRTKCJDQV9kYXRhc2V0Mi5yZHMiKSAj5bey57uP57uP6L+H5YiG57uE5aSE55CG5LqGCmBgYAoKIyMg5L+u5pS55YiG576kCmBgYHtyfQp1bWFwcGxvdChDQV9kYXRhc2V0Miwgc3BsaXQuYnkgPSAic2FtcGxlIikKdW1hcHBsb3QoQ0FfZGF0YXNldDIsZ3JvdXAuYnkgPSAiZ3JvdXBzIiwgc3BsaXQuYnkgPSAiY29uZGl0aW9ucyIpCm11bHRpX2ZlYXR1cmVwbG90KGMoIkhFWTEiLCJHSkE1IiwiU0VNQTNHIiwiQ1hDTDEyIiwiU09YMTciLCJDREg1IiwiUEVDQU0xIiksQ0FfZGF0YXNldDIpCm11bHRpX2ZlYXR1cmVwbG90KGMoIkFDS1IxIiwiUExWQVAiLCJJVEdBNiIsIlBFQ0FNMSIpLENBX2RhdGFzZXQyKQpgYGAKCmBgYHtyfQojIHRhYmxlKENBX2RhdGFzZXQyJHNhbXBsZSkKCiMg5YWz5rOoY2x1c3RlciAzLTYgRUMKIyBjbHVzdGVyIDcg5Z+66LSo57uG6IOeCgptdWx0aV9mZWF0dXJlcGxvdChjKCJMWVoiLCJQVFBSQyIsIkNENjkiLCJFUENBTSIsIkNESDEiLCJQREdGUkIiLCJDT0wxQTIiLCJQRUNBTTEiLCJDTERONSIpLENBX2RhdGFzZXQyKQoKRG90cGxvdChjKCJMVU0iLCJNTVAyIiwiTUdQIiwiRENOIiwiTVlIMTEiLCJBQ1RBMiIsIkNOTjEiLCJUQUdMTiIpLENBX2RhdGFzZXQyKSAjY2x1c3RlciA3IOe7huiDnuWPr+S7peiiq+iupOS4uuaYr21vZHVsYXRlZCBTTUNzCkNBX2RhdGFzZXQyIDwtIEFkZE1vZHVsZVNjb3JlKENBX2RhdGFzZXQyLGxpc3QoYygiTFVNIiwiTU1QMiIsIk1HUCIsIkRDTiIpKSkKQ0FfZGF0YXNldDIgPC0gQWRkTW9kdWxlU2NvcmUoQ0FfZGF0YXNldDIsbGlzdChjKCJNWUgxMSIsIkFDVEEyIiwiQ05OMSIsIlRBR0xOIikpKQptdWx0aV9mZWF0dXJlcGxvdChjKCJDbHVzdGVyMSIsIkxVTSIsIkFDVEEyIiwiVEFHTE4iKSxDQV9kYXRhc2V0MikKbXVsdGlfZmVhdHVyZXBsb3QoYygiTU1QMiIsIkdKQTQiLCJQRUNBTTEiLCJBQ0tSMSIpLCBDQV9kYXRhc2V0MikKbXVsdGlfZmVhdHVyZXBsb3QoYygiTU1QMiIsIkdKQTQiLCJQRUNBTTEiLCJBQ0tSMSIpLCBodW1hbl9jb3JvbmFyeSkKCmBgYAojIOe7huiDnuWIhuexuwpgYGB7cn0KbXVsdGlfZmVhdHVyZXBsb3QoYygiTFlaIiwiUFRQUkMiLCJDRDY5IiwiUERHRlJCIiwiQ09MMUEyIiwiUEVDQU0xIiwiQ0xETjUiKSwgQ0FfZGF0YXNldDEpCmYoIkxVTSIsQ0FfZGF0YXNldDEpCnVtYXBwbG90KENBX2RhdGFzZXQxKQpgYGAKCmBgYHtyfQpDQV9kYXRhc2V0MSAgPC0gQ0FfZGF0YXNldDEgICU+JSBGaW5kTmVpZ2hib3JzKGRpbXMgPSAxOjIwKSAlPiUgRmluZENsdXN0ZXJzKHJlc29sdXRpb24gPSAwLjUpCmBgYAoKCmBgYHtyfQpJZGVudHMoQ0FfZGF0YXNldDEpIDwtIENBX2RhdGFzZXQxJHNldXJhdF9jbHVzdGVycwpDQV9kYXRhc2V0MSA8LSBSZW5hbWVJZGVudHMoQ0FfZGF0YXNldDEsJzAnID0gJ01hY3JvcGhhZ2UnLCcxJz0nU01DJywnMicgPSAnU01DJywnMycgPSAnU01DJywnNCcgPSAnTksgY2VsbCcsCiAgICAgICAgICAgICAgICAgICAgICAgICAgICAnNScgPSAnRUMxJywnNicgPSAnU01DJywnMTAnPSdNYXN0IGNlbGwnLCc4Jz0nTWFjcm9waGFnZScsJzExJz0nTWFjcm9waGFnZScsJzEyJz0nQiBjZWxsJywnNyc9J0VDMScsICc5Jz0nRUMyJykKSWRlbnRzKENBX2RhdGFzZXQxKSA8LSBvcmRlcmVkKElkZW50cyhDQV9kYXRhc2V0MSksbGV2ZWxzID0gYygiQiBjZWxsIiwiRUMxIiwiRUMyIiwiTWFjcm9waGFnZSIsIk1hc3QgY2VsbCIsIk5LIGNlbGwiLCJTTUMiKSkKQ0FfZGF0YXNldDEkQ2xhc3NpZmljYXRpb24xIDwtIElkZW50cyhDQV9kYXRhc2V0MSkKdW1hcHBsb3QoQ0FfZGF0YXNldDEpCiMgc2F2ZVJEUyhDQV9kYXRhc2V0MSwiQ0FfZGF0YXNldDEucmRzIikKCklkZW50cyhDQV9kYXRhc2V0MikgPC0gQ0FfZGF0YXNldDIkc2V1cmF0X2NsdXN0ZXJzCnVtYXBwbG90KENBX2RhdGFzZXQyKQpDQV9kYXRhc2V0MiA8LSBSZW5hbWVJZGVudHMoQ0FfZGF0YXNldDIsJzAnID0gJ1QgY2VsbCcsJzQnID0gJ05LIGNlbGwnLCc5JyA9ICdOZXVyYWwgcHJvZ2VuaXRvcicsJzInID0gJ1N0cm9tYWwgY2VsbCcsJzcnPSdTdHJvbWFsIGNlbGwnLCczJz0nRUMnLCc2Jz0nRUMnLCcxJz0nTWFjcm9waGFnZScsJzUnPSdCIGNlbGwnLCc4Jz0nTWFzdCBjZWxsJywnMTAnPSdCIGNlbGwnKQpJZGVudHMoQ0FfZGF0YXNldDIpIDwtIG9yZGVyZWQoSWRlbnRzKENBX2RhdGFzZXQyKSxsZXZlbHMgPSBjKCJCIGNlbGwiLCdFQycsJ0ZCJywiTWFjcm9waGFnZSIsIk1hc3QgY2VsbCIsJ05ldXJhbCBwcm9nZW5pdG9yJywiTksgY2VsbCIsIlNNQyIsJ1QgY2VsbCcpKQpDQV9kYXRhc2V0MiRDbGFzc2lmaWNhdGlvbjEgPC0gSWRlbnRzKENBX2RhdGFzZXQyKQp1bWFwcGxvdChDQV9kYXRhc2V0MikKCiNzYXZlUkRTKENBX2RhdGFzZXQyLCJDQV9kYXRhc2V0Mi5yZHMiKQoKSWRlbnRzKGh1bWFuX2Nvcm9uYXJ5KSA8LSBodW1hbl9jb3JvbmFyeSRzZXVyYXRfY2x1c3RlcnMKdW1hcHBsb3QoaHVtYW5fY29yb25hcnkpCmh1bWFuX2Nvcm9uYXJ5IDwtIFJlbmFtZUlkZW50cyhodW1hbl9jb3JvbmFyeSwnMCcgPSAnRkInLCcxJyA9ICdNYWNyb3BoYWdlJywnMic9J0VDJywnMyc9J1NNQycsJzQnPSdTTUMnLCc1JyA9ICdUIGNlbGwnLCc2JyA9ICdCIGNlbGwnLCc3Jz0nTmV1cm9uJywnOCc9J1BsYXNtYScsJzknPSdQbGFzbWEnLCcxMCc9J0VDJykKc2F2ZVJEUyhodW1hbl9jb3JvbmFyeSwiaHVtYW5fY29yb25hcnkucmRzIikKQ0FfZGF0YXNldDIkQ2xhc3NpZmljYXRpb24xIDwtIElkZW50cyhDQV9kYXRhc2V0MikKdW1hcHBsb3QoQ0FfZGF0YXNldDIpCmBgYAoKYGBge3J9CkNBX2RhdGFzZXQxX21hcmtlcnMgPC0gRmluZEFsbE1hcmtlcnMoQ0FfZGF0YXNldDEsIGxvZ2ZjLnRocmVzaG9sZCA9IDAuNywgbWluLmRpZmYucGN0ID0gMC4zLG9ubHkucG9zID0gVCkKQ0FfZGF0YXNldDJfbWFya2VycyA8LSBGaW5kQWxsTWFya2VycyhDQV9kYXRhc2V0MiwgbG9nZmMudGhyZXNob2xkID0gMC43LCBtaW4uZGlmZi5wY3QgPSAwLjMsb25seS5wb3MgPSBUKQpodW1hbl9jb3JvbmFyeV9tYXJrZXJzIDwtIEZpbmRBbGxNYXJrZXJzKGh1bWFuX2Nvcm9uYXJ5LCBsb2dmYy50aHJlc2hvbGQgPSAwLjcsIG1pbi5kaWZmLnBjdCA9IDAuMyxvbmx5LnBvcyA9IFQpCmBgYAoKCiMjIOWfuui0qOe7huiDngpgYGB7cn0KZHMwIDwtIHJlYWRSRFMoImRzMC5yZHMiKQpkczEgPC0gcmVhZFJEUygiZHMxLnJkcyIpCmRzMiA8LSByZWFkUkRTKCJkczIucmRzIikKYGBgCgojIyDljrvpmaTovrnop5LnmoTnprvnvqTnu4bog54KYGBge3J9CiMgc2VsZWN0LmNlbGxzIDwtIENlbGxTZWxlY3RvcihwbG90ID0gRGltUGxvdChkczEsIHJlZHVjdGlvbiA9ICJ1bWFwIikpIAojIGRzMSA8LSBzdWJzZXQoZHMxLCBjZWxsID0gc2VsZWN0LmNlbGxzKQojIHVtYXBwbG90KGRzMSkKIyBzYXZlUkRTKGRzMSwiZHMxLnJkcyIpCgoKIyBzZWxlY3QuY2VsbHMgPC0gQ2VsbFNlbGVjdG9yKHBsb3QgPSBEaW1QbG90KGRzMCwgcmVkdWN0aW9uID0gInVtYXAiKSkgCiMgZHMwIDwtIHN1YnNldChkczAsIGNlbGwgPSBzZWxlY3QuY2VsbHMpCiMgdW1hcHBsb3QoZHMwKQojIHNhdmVSRFMoZHMwLCJkczAucmRzIikKYGBgCgojIyDnu4bliIbkuprnvqQKYGBge3J9CklkZW50cyhkczIpIDwtIGRzMiRjb25kaXRpb25zCmRzMl9BQyA8LSBzdWJzZXQoZHMyLCBpZGVudHMgPSAiQUMiKQpkczJfUEEgPC0gc3Vic2V0KGRzMiwgaWRlbnRzID0gIlBBIikKCmRzMCAgPC0gZHMwICU+JSBGaW5kTmVpZ2hib3JzKGRpbXMgPSAxOjIwKSAlPiUgRmluZENsdXN0ZXJzKHJlc29sdXRpb24gPSAwLjIpCmRzMSAgPC0gZHMxICU+JSBGaW5kTmVpZ2hib3JzKGRpbXMgPSAxOjIwKSAlPiUgRmluZENsdXN0ZXJzKHJlc29sdXRpb24gPSAwLjIpCgpkczJfQUMgPC0gZHMyX0FDICU+JSBGaW5kTmVpZ2hib3JzKGRpbXMgPSAxOjIwKSAlPiUgRmluZENsdXN0ZXJzKHJlc29sdXRpb24gPSAwLjE1KQpkczJfUEEgPC0gZHMyX1BBICU+JSBGaW5kTmVpZ2hib3JzKGRpbXMgPSAxOjIwKSAlPiUgRmluZENsdXN0ZXJzKHJlc29sdXRpb24gPSAwLjEpCgpkczJfQUNfbWFya2VycyA8LSBGaW5kQWxsTWFya2VycyhkczJfQUMsIGxvZ2ZjLnRocmVzaG9sZCA9IDAuNywgbWluLmRpZmYucGN0ID0gMC4zLG9ubHkucG9zID0gVCkKZHMyX1BBX21hcmtlcnMgPC0gRmluZEFsbE1hcmtlcnMoZHMyX1BBLCBsb2dmYy50aHJlc2hvbGQgPSAwLjcsIG1pbi5kaWZmLnBjdCA9IDAuMyxvbmx5LnBvcyA9IFQpCmRzMG1hcmtlcnMgPC0gRmluZEFsbE1hcmtlcnMoZHMwLCBsb2dmYy50aHJlc2hvbGQgPSAwLjcsIG1pbi5kaWZmLnBjdCA9IDAuMyxvbmx5LnBvcyA9IFQpCmRzMW1hcmtlcnMgPC0gRmluZEFsbE1hcmtlcnMoZHMxLCBsb2dmYy50aHJlc2hvbGQgPSAwLjcsIG1pbi5kaWZmLnBjdCA9IDAuMyxvbmx5LnBvcyA9IFQpCmRzMm1hcmtlcnMgPC0gRmluZEFsbE1hcmtlcnMoZHMyLCBsb2dmYy50aHJlc2hvbGQgPSAwLjcsIG1pbi5kaWZmLnBjdCA9IDAuMyxvbmx5LnBvcyA9IFQpCgp1bWFwcGxvdChkczApCnVtYXBwbG90KGRzMSkKCnVtYXBwbG90KGRzMl9BQykKdW1hcHBsb3QoZHMyX1BBKQojdW1hcHBsb3QoZHMyLHNwbGl0LmJ5ID0gImNvbmRpdGlvbnMiLGdyb3VwLmJ5ID0gIkNsYXNzaWZpY2F0aW9uMSIpCgoKSWRlbnRzKGRzMCkgPC0gZHMwJHNldXJhdF9jbHVzdGVycwpkczAgPC0gUmVuYW1lSWRlbnRzKGRzMCwnMCcgPSAnRmlicm9ibGFzdCcsJzEnPSdTTUMnLCcyJyA9ICdGaWJyb215b2N5dGUnLCczJyA9ICdQZXJpY3l0ZScsJzQnPSdQZXJpY3l0ZScsJzUnPSdVbmtub3duJykKCnVtYXBwbG90KGRzMCkKCklkZW50cyhkczEpIDwtIGRzMSRzZXVyYXRfY2x1c3RlcnMKZHMxIDwtIFJlbmFtZUlkZW50cyhkczEsJzAnID0gJ1NNQzEnLCcxJz0nRmlicm9teW9jeXRlJywnMicgPSAnVW5rb3duMScsJzMnID0gJ1NNQzInLCc0Jz0nVW5rbm93bjInKQoKSWRlbnRzKGRzMikgPC0gZHMyJHNldXJhdF9jbHVzdGVycwpkczIgPC0gUmVuYW1lSWRlbnRzKGRzMiwnMCcgPSAnU01DMScsJzEnPSdGaWJyb215b2N5dGUnLCcyJyA9ICdTTUMyJywnMycgPSAnRkInLCc0Jz0nU01DMycpCgp1bWFwcGxvdChkczIpCnVtYXBwbG90KGRzMSkKYGBgCgpgYGB7cn0KSWRlbnRzKGRzMl9QQSkgPC0gZHMyX1BBJHNldXJhdF9jbHVzdGVycwp1bWFwcGxvdChkczJfUEEpCmRzMl9QQSA8LSBSZW5hbWVJZGVudHMoZHMyX1BBLCcwJz0nRmlicm9teW9jeXRlJywnMScgPSAnU01DMScsJzInID0gJ1NNQzInKQpzYXZlUkRTKGRzMl9QQSwgZmlsZSA9ICJkczJfUEEucmRzIikKCklkZW50cyhkczJfQUMpIDwtIGRzMl9BQyRzZXVyYXRfY2x1c3RlcnMKdW1hcHBsb3QoZHMyX0FDKQpkczJfQUMgPC0gUmVuYW1lSWRlbnRzKGRzMl9BQywnMCc9J0ZpYnJvYmxhc3QnLCcxJyA9ICdTTUMxJywnMycgPSAnU01DMicsJzInID0gJ0ZpYnJvbXlvY3l0ZScpCklkZW50cyhkczJfQUMpIDwtIGZhY3RvcihJZGVudHMoZHMyX0FDKSxsZXZlbHMgPSBjKCJGaWJyb215b2N5dGUiLCJTTUMxIiwiU01DMiIsIkZpYnJvYmxhc3QiKSkKc2F2ZVJEUyhkczJfQUMsIGZpbGUgPSAiZHMyX0FDLnJkcyIpCgpnZ3NhdmUoImRzMl9QQXVtYXAuc3ZnIixwbG90ID0gdW1hcHBsb3QoZHMyX1BBKSxkZXZpY2UgPSBzdmcsIHdpZHRoID0gNiwgaGVpZ2h0ID0gNCkKZ2dzYXZlKCJkczJfQUN1bWFwLnN2ZyIscGxvdCA9IHVtYXBwbG90KGRzMl9BQyksZGV2aWNlID0gc3ZnLCB3aWR0aCA9IDYsIGhlaWdodCA9IDQpCmBgYAoKYGBge3J9CndyaXRlLmNzdihkczJfQUNfbWFya2VycywiZHMyX0FDX21hcmtlcnMuY3N2IikKd3JpdGUuY3N2KGRzMl9QQV9tYXJrZXJzLCJkczJfUEFfbWFya2Vycy5jc3YiKQoKYGBgCgo=